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1 Introduction 

In recent years much attention lias been devoted to tlie study of mean field 
equations of Liouville type on Riemann surfaces and in the presence of singular 
sources. Such an interest has originated from various areas of mathematics 
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and physics, starting with the assigned Gauss curvature problem which can 
be reduced to analyze the solution set of 

Au + K{x) e^" = in (1) 



for a given function K defined in R^. If m is a solution to ([T]), then the metric 
g = e^** is conformal to the flat metric |(ia;p and such that K is the Gaus- 
sian curvature of the new metric g. Equation ([1]) also appears in the analysis 
of gravitating systems, in the statistical mechanics description of the vorticity 
in fluid mechanics (see [2|3Pf20l[2T] ) and has been studied more recently in 
the context of self-dual gauge field vortices (see [TBll35f32j ). 

The solution set of equation ([T]) depends very much on the properties of K. 
When K is negative, uniqueness results are always available, while for total 
positive curvature K, either uniqueness or multiplicity of solutions holds, de- 
pending on K. See [2lJ for various examples. In this paper we will focus on 
particular cases of positive functions K. 

We are also interested in the problem of the assigned Gauss curvature with 
conical singularities (see [5Hl33|27|14|lUj ). For a given Riemann surface (M, g), 
we aim at determining the range of the parameters A, p G M such that 



is solvable on M, where is the Laplace-Beltrami operator, dag is the volume 
element corresponding to the metric g, f E C{M) with Jm f dcg = and 6p 
is the Dirac measure with singularity at P G M. In case M has a non-empty 
boundary, both Dirichlet or Neumann boundary conditions on dM are of 
interest for the applications. One could also consider sums of Dirac measures 
located at several source points. 

In the applications, the simplest situations correspond to the 2-sphere M = 
and the fiat 2-torus M = C/ (^iZ + ^2^), with periodic cell domain generated 
by (^1 and ^2- Recall that the 2-sphere, with the standard metric induced by 
Lebesgue's measure in R^, has already played a special role in the assigned 
curvature problem (see [T7|19|18j ). while the torus is important since many 
vortex-like configurations naturally develop into periodic lattices. 

For the sphere, p = corresponds to a particular case of the so-called Onsager 
vortex problem (see [MIES])- In fact, for closed surfaces and when there is 
no singularity {i.e. dM = and p = 0), the solvability of ([2]) is quite well 
understood in terms of the topological properties of M. Starting with the 
work of Y. Y. Li in [23], subsequently completed by C.-C. Chen and C.-S. 
Lin in [51IS] . we know that, when p = 0, the solutions of ([2]) are uniformly 



2 



bounded for any fixed A G M \ 47rN, and the Leray-Schauder degree dx of 
the corresponding Fredholm operator can be exphcitly computed: for A G 
(47r(m - l),47rm), m G N*, 

, . .f . , (-x(M) + l) ■■■ (-;^(M)+m-l) 

cJa = 1 if m = 1 , dx = 7 -T^ if m > 2 , 

(m — 1)! 



where x{^) is the Euler characteristics of M (see p[25] ). Actually, for the flat 
2-torus we have that dx = 1 also when A G 47rN, and so ([2]) has a solution for 
every A G M, if p = (see [6j). For the standard 2-sphere, we have that dx = ^ 
for all A > Stt. But, by a more precise topological argument (see [13]), it can 
be shown that in this case equation ([2]) admits a solution for any A G M \ 47rN, 
if p = 0, and also some multiplicity results can be proved. 

The situation is much more complex in presence of a Dirac measure. In fact, 
an expression for the Leray-Schauder degree, given by dx = — x(^) + 2, is 
available only when p > 1 and A G (47r, Svr) (see [7]). Hence, for S^, such a 
degree formula yields no information about the solvability of ([2]), and indeed 
this issue turns out to be very delicate (see [Mf33] ). A similarly delicate situ- 
ation occurs for the flat 2-torus when A = 47r and p = 2. Consider for instance 
the equation 

AgU + ATiiy^ ^J^^da -^p)=^ inM = C/(^iZ + 6Z), 

with P = 0. Then, C.-S. Lin and C. L. Wang have shown in p6j that there is 
no solution in the case of a rectangular lattice {i.e. = a, C,2 = ib, a, b > 0), 
while there is a solution for a rhombus lattice {i.e. = a, ^2 = ae*'^/^, a > 0). 

As usual, the cause for such surprising existence or nonexistence situations, 
there is a lack of compactness for the solution set. It can explained by looking 
at the singular Liouville equation 

Aw + e^" = 27rp5^o inR^ (3) 



The solutions of ([3]) are described using complex notations (with R^ ~ C), by 
means of the Liouville formula: 



u{z) 




4|f(^)P ' 
2^2 



1 + 



where / is a meromorphic function. To account for the singularity at 2:05 "we 
must require that it corresponds to a pole of /, or a zero of /', or a branch 
point of /, with the appropriate order. In particular, for all solutions u of 
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that satisfy e^" G L^(R2), there holds: 

J e=^" dx = 47r (1 + p) , 

(see [8]|9|29] ). This explains the role of Atc for the regular problem, p = 0. 

One of the motivations of this paper is to connect the solvability of to some 
weighted Liouville-type equations in which generalize ([3]). To see how this 
class of equations could arise, let us first focus on problem ([2]) over the square 
2-torus. In this case we know that a solution exists if we replace the Dirac 
measure by a smooth function (see [6]). Let Q = (—1,1)^ C and suppose 
that the source point P coincides with the origin G For given A > and 
p > — 1, denote by the solution of 

j + ' [j^. - m) ° ( m^'' ' W\) (4) 

I Ue doubly periodic on dil , 



where Cg := Jq ^^^2^|^|2-)2 dx converges to 1 as e ^ 0. By doubly periodic on 
dQ, we mean u{x + 2e) = u{x) for any x G dQ, with e = (1, 0) or e = (0, 1). 
To prove the existence of a solution of ([2]), it is natural to investigate under 
which conditions on A and p we can pass to the limit in (jl]), along with a 
subsequence of £ ^ 0. In other words, whenever possible, we need to establish 
a priori estimates for in suitable norms. To this end, denote by m^^q the 
unique solution to the problem 

[ u doubly periodic on dil , JqU dx = , 

which takes the form Ueflix) = | log(£^ + \x\'^) + ipei^), for some suitable 
function ip^, which is uniformly bounded in C^''^-norm, with respect to e > 0. 
Then = + Uefl is a solution of (jl]) if and only if satisfies: 

e2(«,.o+..) dx ~ My " ° ^ ' 

f £ doubly periodic on dVL . 



Av, + A 



The function e'^^ " is bounded from above and from below away from zero in 
C^Q^{Q \ {0}), uniformly in e G (0,1). Therefore, by well known estimates 
based on blow-up analysis (see, e.g., [23ll5]|6][T|3T] ). we know that for A ^ 47rN, 
then Ue is bounded uniformly in C?'Z{^ \ {0}) with respect to £ > 0. 
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Next, we have to investigate what may happen in a neighborhood of the origin 
G fi. To this purpose, define := + | log A — | log {J^ e^^^ dx) and 
:= e^^^^ in -8^(0) C Q, for r > small enough. There holds: 



J ^ ' ' ^ ^ ^ Le^^'^dx ~ 

Br{0) 



If We was not uniformly bounded in Br{0), then we could find a sequence 
(£n)neN with lim„^oo£„ = 0, and a sequence of points in 3^(0) with 



lim„_>oo Xn — such that 



We (Xn) — maXU^e 



+00 as n — > oo , 



(5) 



For any n E N, define s„ := max{e„, | x„ | , exp ( — ^2(i+p) ) } ^^d observe that 
jirn^Sn = 0. Let Rn := We„{xn + Snx), Bn := Sr/s„(0). For n large, C/n(x) := 
WsrX^n + Snx) + 2{1 + p) log s„ satisfics 



+ 



2^ + a; 



Rne^^- + o{l) in Br,, 



Un{0) = + 2 (1 + p) logs„ , 



+ 



Xr. 



+ X 



2\ P 



Rn e^^" dx<X. 



By definition of we know that limsup^^^ s„ exp ^^^^^^^^^ ) > 1. We do not 
know whether this limit is finite or not. If 



1- / WeS^n) \ ^ 

hmsup Sn exp — — — - < oo , 



(6) 



using Harnack's estimates, we can determine a subsequence along which 

Un ^ c/oo in q'oc(^') ' - ^ £oo e [0, 1] , — ^ e s(o, i) c 
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and Ws„{xn + pointwise converges to a positive constant Woo > 0. In 
addition, Uoo satisfies: 

-Af/oo = {el + \xoo + x\y Woo e^^- in ^ 
< f/oo(0) = max t/oo > , 
I + |a;oo + x\y Woo e'^°° < A . 

If £00 = 0, then by well known classification results (see p[29] ). we know that 
Woo J \xoo + xl'^f e^^=° dx = 47r (1 + p) . 

R2 

So, we could rule out the occurrence of ([5]) in this case by restricting the 
problem to A < 47r (1 + p). 

Hence, assume that £00 > 0. Then by adding a constant to U, scaling and 
translating we arrive at the limiting problem: 

-Af/ = (1 + IxH^e^^ inM^^ (a) 
' l{l + \x\ye^^dx<\. (b) ^'^^ 

.R2 



Now, with £00 > 0, it is much harder to identify the range of values as- 
sumed by /]k2(1 + \x\'^)'' e^'^dx for all solutions of ([7^). It is no longer a single- 
ton, but a whole interval contained in {2tt (1 + p),4:'n' (1 + p)) and containing 
(47r max(l, p), 47r (1 + p)), whose explicit range is still under investigation, 
even when we restrict the problem to radially symmetric solutions of ([7^), see 
[TT|24|25j . In this paper, we are going to identify necessary conditions on A 
that allow us to rule out the occurrence of ([5]) under Condition Qj. In other 
words, if, for some Aq, there is no solution of ([7]) for any A < Aq, then ([2]) with 
M = Q would have a solution, provided Condition holds. 

On the other hand, consider on the solutions of 

-Am = (1 + IxH^e^" inR2, 
' J{l + \x\Ye'''dx = X. 

,R2 
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This problem is equivalent to ([2]). Let S : — be the stereographic 
projection with respect to the north pole, N := (0,0,1). To any solution u 
of dHD, we can associate 

v{y) ■.= u{x) -l-\og{ / )-log2, (9) 



1+ X 



for X = Tj(y) and any y G § . The function v solves 



Moreover, if m is a locally bounded solution of ([8]), according to [llj, u{x) 
— ^ log |x| as |x| oo, which shows that 



y&2,y-,N \og\y — N\ 2n 



Here, by |?/ — N|, we denote the euclidean distance from ?/ to N in D S^. As 
a consequence. 



Agv = (a - 27r (iV + 2)) 6m + 0(1) 



as y ^ M. Hence w is a solution of ([2]) in S"^ with p = — {N + 2) and / = 0. 
Here the parameter p is not the same as in ([7]). Notice that v given by (Q is 
such that /g2 e'^'" dag = 1. As a special case, if A = 27r (A^ + 2), we find that v is 
a bounded solution of ([2]) with p = 0. Viceversa, if f is a solution of ([2]) with 
/ = 0, then 

M = i; o S"' - - log I / e''' dag 1 + log 2 + - log ^ 




2 ^ \ (1 + |x|2)^+2 



is a solution of ([HD- So, existence and multiplicity of solutions to ([HD provide 
us with existence and multiplicity of solutions to ([2]) on the 2-sphere. 

From now on, we shall focus on the study of ([8]), for all positive values of A^. 
Our main result establishes the existence of many bounded, radial solutions. 

Theorem 1 For allk>2 and N > k{k + l)-2, there are at least 2{k-2) + 2 
distinct radial solutions of ^ with A = 2ti {N + 2), one of them being the 
function ulj{r) := i log L^l^t^U ) with r =\x\. 
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To conclude this introduction let us mention that in many papers devoted 
to other applications than assigned Gauss curvature problems, the conven- 
tion is to consider equation ([2]) with e" instead of e^". Necessary adaptations 
are straightforward (replace -u by | (m — log 2)) and therefore are left to the 
interested reader. 



2 Known results 



In this section, we collect known existence results and some of the properties 
of the bounded solutions to 



Am + (1 + |xP)V" = 0, xe 



and + IxH^^rfx < oo, (11) 



where iV > is a real parameter which in the applications usually enters as 
an integer value. Several papers, mostly by C.-S. Lin et al. |TT|24f25j . have 
dealt with this class of equations. Almost all results are concerned with the 
set of radially symmetric solutions of fill I) , which can be parametrized by a 
parameter a G R as follows: 

< + V + (1 + r2)^e2"- = in (0, +oo) , 

7 (12) 
Ma(0)=a, m;^(0)=0, [l + r"^)^ e^'''^ r dr < +O0 . 



It is proven in [TT|24ll25j that, for every a G M, there exists a unique solution Ua 
of ffT2l). which moreover satisfies 



\mv^{ua{r) + a{a) logr) = (3{a) , (13) 
where a ^— a(a) and a ^ (3{a) are two C^(M, M) functions defined by: 

oo oo 

y(l + r2)^e^"-rdr and (3{a) = J {1 + r^)^ e^""" r logr dr . 



oo 

a(a] 



Moreover it was proved in [TT] (also see [25]) that, for any > 0, 

lim a{a) = 2{N + 1) and lim a{a) = 2 minil, N} . (14) 

a — ^— oo a — *-+cxD 
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Pohozaev's identity applied to (fTTj) shows that a{a) G {2,2{N + 1)). For 
integrabihty reasons, we also know that a{a) > + 1, and so 



max{2,iV+ 1} < a(a) < 2(A^+ 1) VaGM. 

Our problem is to find the solutions of (fT2|) corresponding to a given a, that is 
u" + ^ + {l + r2)^e2" = in (0, +oo) , 
m'(0) = 0, [{l + rYe^''rdr = a. 







In [2H[25] . C.-S. Lin investigates the uniqueness issue for problem (fT5l) by 
identifying the values of a for which there is a unique a such that a{a) = a. 
He proves in [25j that uniqueness holds for a G (2A^, 2(A^ + 1)) if > 1 and 
for all a G (2, 2(A^ + 1)) if iV < 1. On the other hand, it is easy to verify that 
for all A^, the function 



is a solution to ( |T2|) for := I log (2(A^ + 2)), and satisfies a{a*^) = 

N + 2. Since + 2 < 2A^ < 2(A^ + 1) for all A^ > 2, by continuity of 
a ^ a{a) , it appears that there exists at least two different values of a such 
that a{a) = a, for any a G (miUagiR a{a), 2N). In other words, for those values 
of a there exists at least two radially symmetric solutions of (fTTj) satisfying: 

j{l + \x\Ye^" dx = 27ra. 

R2 



Moreover for N = 2, 0(02) = 4 = 2A^ and a'(ci2) < 0, so the above multiplicity 
results also holds true for N = 2. Summarizing, we can state the following 
result. 

Theorem 2 p!T|24|25j Let N be any positive real number. 

(i) If N < 1, then the curve a t— > a{a) is monotone decreasing. Moreover, 
there exists a radially symmetric solution u of / f73]) if and only if a E 
(2,2(A^+ 1)), and such a solution is unique. 

(ii) If N > 1, then for all a G {2N,2{N + 1)), there exists a unique a G M 
such that a{a) = a. In other words, for such a, problem [TS\) is satisfied 
by a unique radial solution. 
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(iii) If N > 2, then minagK a(a) < 2N, and for all a G (miiiagK a;(a), 2N), 
there exists at least two radial solutions of l[TS\) . 

Remark. Concerning part (iii) of Theorem [21 by a closer inspection of the 
results of [21], actually we know that for N > 2, problem f[T^ is satisfied by 
a unique radial solution also when a = 2N. 

An important tool in the proof of the above uniqueness results is the study of 
the linearized problem 

f + ^ + 2 (1 + rXe^""^ <^a = , r e (0, +oo) , 

(17) 

U,(0) = 1, v>'aiO) = 0, 



and in particular the number of zeroes of the function ipa when a is in the range 
{2N,2{N + 1)). The number of critical points of a a{a) is also connected 
with the number of zeroes of ipa in the range (minagKa(a), 2A^). It is indeed 
easy to prove that as r goes to +oo, 

ipa{r) ~ -a'{a) logr + b'{a) + o(l) , (18) 



and hence, that <^a is a bounded function if and only if a G M is a critical 
point of the function a. As a special case, for all A^, if miUagK a{a) is achieved 
for some finite a, then ipa is bounded. 



2.1 Non radially symmetric solutions 

By (El), to any solution Ua of (lT2l) . we can associate a function Va on S^, such 
that /g2 e'^'" dag = 1, which solves (ITOl) for A = 27r a;(a). At level a = N + 2, Va 
is a bounded solution of (121) (with / = 0, p = 0), which is axially symmetric 
with respect to the unit vector (0, 0, 1) pointing towards the north pole N 
of S^. Since = Va*^ is the unique constant solution of ([21), if we know the 
existence of more than one solution at level a = A^ + 2, then there is an axially 
symmetric solution of f[TUl) which is not constant, and that can be thus rotated 
in order to be axially symmetric with respect to any vector e G \ {N, S}. 
Let us denote by such a solution. Applying ([9D to v^, we find a solution 
of (ITTl) which is not radially symmetric. If at level A = 2tt{N + 2) we find k 
solutions of (ITTl) different from m^, then we get k punctured spheres of non 
radially symmetric solutions of (ITTl) . Details can be found in [25] . 
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3 Multiplicity results for radially symmetric solutions 

For given A^, critical levels of the curve a ^ a{a) determine the multiplicity 
of the radial solutions at a given level. The number of zeros of the solutions 
of the linearized problem can change only at critical points of a, see below 
Section [3^ In the special case a = + 2, a bifurcation argument provides us 
with a very precise multiplicity result, which is our main result, see Section l373l 

3.1 A preliminary result 

Let {c^)lz\ be the ordered sequence of all critical values, counted with mul- 
tiplicity, of the curve a a{a), Cq := infagjR a (a) and = 2(A^ + 1). Denote 
by {a^)'^z\ a sequence of critical points corresponding to {c^)'^z\ and, for any 
= 1, 2, . . . n — 1, let := +2 if is a local minimum, := —2 if is a 
local maximum, and := otherwise. Also let := 2 if infaGRa(a) is not 
achieved, and otherwise. Let Xn{o) := 1 if a > 2A^ and otherwise. The 
next proposition links these values with the number of solutions of f|T5l) . 

Proposition 3 Let N he any positive real number. With the above notations, 
for any a > 0, Equation 173]) has exactly Z]f=o ~ xi^^) solutions such that 
a{a) = aifae (c^,cf+i) n (M\ {2A^}), for any k = 0, l,...n-l. 

The proof is straightforward and left to the reader. We shall now focus on 
the study of the critical points of the curve a a{a). Our main tool is the 
linearization of f|T2|) . 

3.2 Study of the linearized problem 

In order to study the multiplicity of radial solutions for ([H]), it is convenient 
to perform the Emden- Fowler transformation in the linearized equation (ITTl) : 

t = logr , Wa{t) := ipa{r) ■ 

The equation in f|T7j) is then transformed into 

+ 2 e^\l + e2*)^e2""(^') Wait) = , tE (-oo, +oo) . (19) 

When the equation for w% := Wa* reads 

^rW + ^^^^^^W = 0, tG(-oo,+oo). (20) 
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With one more change of variables, w{t) = ip{s), s = tanht, we find Legendre's 
equation: 

A (^(1 _ ^\ ^ ^ + 2 ^ ^ Q 

ds \ ds j 2 

which defines the Legendre polynomial of order A; G N* if + 2 = k{k + 1). 
Notice that the composition of the two above changes of variables amounts to 
write s = ^J^qr^' ^ very precise spectral analysis made in [22] shows that the 
above equation has bounded solutions if and only if there is a positive integer 
k such that 1 + 2A; = y^l + 4(A^ + 2), that is, if and only if 

^ ^ -1 + 7l + 4(iV + 2) 
9T(iV) := ^ ^ 

is a positive integer. This is solved by = A^^ := k{k + 1) — 2, A; G N*. 
Actually, we are interested only in /c > 2 since we only deal with A^ > 0. As 
a consequence, is a critical point of a if and only if Dl(A^) is a positive 
integer: A^ = 4, 10, 18, . . . 

For A^ = A^fc > 0, we know explicitly the solutions to ( l20l) . They are the 
Legendre polynomials, namely, with s = tanht, 

w%^{t) = Pk{s) for all integer k>2. (21) 

Lemma 4 [22] Take N > 1. Then, there exist hounded solutions of [2D\) if 
and only if DT(A^) is a positive integer. In such a case, ipa*^ has exactly ^{N) 
zeroes in the interval (— oo, +oo). 

For all a G M, has at least two zeroes in the interval (0, +cxd). This ob- 
servation is a key step in the uniqueness proofs of [Mf25] . Zeroes of (fa will 
also play an important role in multiplicity results. The next observation is a 
standard result for linear ordinary differential equations. 

Lemma 5 For any A^ > 0, Oq > and R > 0, if ipao has k zeroes in (0, R) 
and ifa^iR) 7^ 0, then there exists an e > such that (fa o-lso has exactly k 
zeroes in (0, R) for any a G (ao — aQ + e). 

Proof. If the result were wrong, then we could find some r G (0, i?) such that 
(faoir) = f'aoir) = and so we would get v^ao = 0. □ 

Corollary 6 For any N G [Nk, Nk+i) , k E N, k > 1, solutions to [2^1 have 
exactly k + 1 zeroes in the interval (—00, +cx)). 
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Proof. If we normalize all functions w^j so that \imt^^oo{W^{t),w'^' (t)) = 
(1, 0), then by continuity in the number of zeros changes if and only if 
is bounded, i.e. if DT(A^) is a positive integer, by Lemma HI At A^i = 0, we 
have: w^lt) = tanh(t). A careful analysis shows that, as a function of A^, the 
number of zeroes is continuous from the right and increasing. □ 

Lemma 7 Take N > 1 and consider ai, 02 G M such that a'{ai) = a'{ai) = 
and a'{a) 7^ if a e (01,02). Then, for all a e (01,02), the functions have 
the same number of zeroes. 

Proof. The proof is based on the same arguments as the proof of Corollary [61 
One has just to replace the continuity in A^ by the continuity in o. □ 

Lemma E] means that when o varies, zeroes may appear or disappear only at 
infinity. For a given o, the sign of the function 

+00 

Jjv(a) :=/(! + rYe^'^^'vl r dr (22) 


governs the dynamics of the zeroes of at infinity as follows. Denote by 
r(o) := max{r > : fair) = 0} the largest zero of (pa- 

Lemma 8 Let a > be such that, for ( = ±1, liTaa-,a,({a-a)>oi^{0') = 
Then there exists e > such that, on (a — e,d) if ( = —1, on (o, a + e) if 

C = +i, 

r(a) 

— {a) = - j {l + r'^ f e^^'^if^rdr 

da r(o) |(^'„(r(o))P 7 

and ^ Jn{o) < if JN{d) ^ 0. 

Proof. First, we choose e > Q small enough so that/o^"^(l+r2)^e2«'" ipl r dr and 
Jn{(i) have the same sign, if Jn{(i) 7^ 0. Next, we take 6 > 0, small. Multiplying 
the equation satisfied by ipa and ipa+b by r (fa+b and r ipa respectively, and 
integrating by parts in the interval (0,r(o)), we get 

r(a) 

2 J {l + rY{e^^'^-e^^^+'')fa¥>a+brdr = -r{a)f',{r{a))y^a+b{r{a)). 


By definition of (fa and using the uniform continuity properties of the functions 
Ua and ipa on (0,r(o)), we obtain: \\ua+b-Ua - bipa\\L^io,ria)) = o{h) as 6 ^ 0, 
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(1 _ e^''^^) ~ -2b<fa- Since <fa+bir{a)) = <f'^{r{a)){r{a + b)) - r{a)) + o{b), we 
get, as 6 ^ 0, 



r(a) 



2 r{a + b)) — r(a) 




•lrdr = r(a) |^l(r(a))| 



6 



+ 0(6) . 







□ 



Notice that, with the notations of Lemma [HI <^a is a bounded function. As a 
consequence, we have the following result. 

Corollary 9 Let a be a critical point of a. There exists e > 0, small enough, 
such that the following properties hold. 

(i) If Jn{cl) > and if, for any a E {a — e, d), all functions (pa are unbounded 
and have k zeroes in (0, +cxo), then ipa is bounded and has k zeroes, and 
for any a E {a, a + e) , (pa is unbounded and has either k or k + 1 zeroes 
in (0, +oo). 

(ii) If Jiy{d) < and if, for any a G {a — e, d), all functions ipa are unbounded 
and have k zeroes in (0, +C)o), then (fa is bounded and has either k ork — 1 
zeroes, and for any a G (a, d+e), ipa has the same number of zeroes as ipa- 

As already seen in Proposition [3l if a is a local extremum of a, the number 
of zeroes changes when a goes through a, since the sign of ipa at infinity also 
changes, by (ITSll . Otherwise, if d is an inflection point, the number of zeroes is 
constant when a passes through d. This explains the ambiguity in the previous 
result. 

Actually for the particular N = Nk we can exactly compute the 

value of Jn as follows: 

Proposition 10 Let us define j{k) := Jn{(i*nJ for any integer k > 2. Then, 
j{k) = if k is odd, and j{k) > if k even. 

Proof. By using (1^ . we can easily compute 



Now, when k is odd, is also odd and so, j{k) = 0. On the contrary, Gaunt's 
formula, see Identity (14), page 195] shows that j{k) > if A; is even. □ 




(23) 



-1 
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Remark. If = Nk, k even, is a local minimum for the fonction a. 
Indeed, this follows from the fact that we know that JAr(a^) > and the 
number of zeroes of Pk- From this we infer that the function is positive at 
infinity for N close to with N < N^, and negative at infinity for close 
to Nk with N > Nk- Hence, a is decreasing to the left of and increasing 
afterwards. 



3.3 A multiplicity result at level a = N + 2 

This section is devoted to the proof of Theorem [H that is a multiplicity result 
for the solutions of problem ( JTSl) at level a = N + 2, which also helps to illus- 
trate Theorem 1.10 in [llj. As seen in the introduction this amounts to study 
the number of bounded solutions to ([2]) with p = (without singularities). We 
show that when X = 2tt a and a = A^ + 2 becomes large, there are more and 
more bounded solutions to ([2]). 

Let M be a radial solution of ([8]) with a = A^ + 2. We may reformulate this 
problem in terms oi f := u — u*^ G P^'^(R^) as a solution to: 

VII/ ^ I I / 

with ^ = 2(A^ + 2). Solutions of ([24]) are bounded by Moreover, ([24]) is 
trivially invariant under the Kelvin transformation: 

Lemma 11 If f is a solution of [24\ ), then the function x y-^ f is also 

a solution of (24^ . 

This lemma allows us to characterize many branches of solutions of (12^ . 

Theorem 12 The function f = is a trivial solution of [24\ ) for any fi > 0. 
For any k >2, there are two continuous half-branches, andC^ , of solutions 
{fi, f) of [24\ ) bifurcating from the branch of trivial solutions at, and only at, 
(/ifc = 2k(k + 1), 0). Solutions in are such that ±/(0) > 0. 

Away from the trivial solutions, all branches are disjoint, unbounded and char- 
acterized by the number of zeroes. In C^, the solutions of ([^]) have exactly k 
zeroes. If k is odd, the branch is the image of by the Kelvin transform. 
If k is even, the half-branches are invariant under the Kelvin transform. 
Finally, for k > 3 and are locally bounded in /i. 

We divide the proof in three steps. 
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Step 1: Existence of unbounded branches of solutions. 

We use a bifurcation method to study the set of radial solutions of (1241) with 
the bifurcation parameter /i. Since branches may be multi-valued in terms 
of A^, we will reparametrize them with a parameter s. Classical results apply 
for instance in R+ x (P^-^ n C^), see p^llHO] . 

By Lemma m there is local bifurcation from the trivial line {(/x, 0)} at the 
points (/ifc := 2{Nk + 2),0), and there is no other bifurcation point in this 
trivial branch. By the properties of the Legendre polynomials P^, if we denote 
by the two continuous half-branches of non-trivial solutions that meet at 
(/ifc, 0), one easily proves that for any (/x, /) G in a neighborhood of (/ifc, 0), 
/(r) has exactly k zeroes in the interval (0, -|-oo). 

Actually, on C^, the number of zeroes of the solutions is constant, namely 
equal to k. For instance, let us prove it for by smoothly parametrizing the 
branch as follows: 

with /i(0) = 2(iVfc + 2) , and / = /, with /o = , sgM. (25) 

Let Afc:={sG(0,oo) : fs admits exactly k zeroes}. Clearly is not empty. 
Since a solution / ^ of fl2^ cannot vanish at a point together with its 
derivative, the smoothness of the map s ^ fs ensures that A^ is open. We 
check that A^ is also closed. Indeed, for ^ A^ with lim^^oo Sn = E 

(0, oo), we see that the zeroes of = fs„ cannot collide at a point vq G M"*" 
since there, lim„^oo fn = fsa^ ^ would vanish together with its derivatives. 
So, for a zero to appear or disappear, this would require the existence of a 
sequence (r„)neN with fnirn) = and lim„^oo''"n = oo. But then, /(r) = 
/s_^(l/r), which is still a solution of fl24p by Lemma [TTl would vanish at 
r = together with its derivatives, and this is again impossible. In conclusion, 
Afc = (0,+cx)). 

Consequently, non-trivial branches with different k cannot intersect or join 
two different points of bifurcation in the trivial branch. For any k > 2, the 
half-branches are therefore unbounded and we can distinguish them as 
follows: if (/X,/) G C^, resp. (/i, /) G C^, then /(O) > 0, resp. /(O) < 0. 

Step 2: Symmetry under Kelvin transform 

Branches of solutions of fl2^ have an interesting symmetry property. If k is 
odd, the solutions in the branches have an odd number of zeroes in (0, +oo) 
and so, they cannot be invariant under the Kelvin transform, because they take 
values of different sign at and near +oo. Since /i^ = 2{Nk -|- 2) is a simple 
bifurcation point, the only possibility is that the branches transform into 
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each other through the Kelvin transform. Otherwise, there would be at least 
four half-branches bifurcating from (yUfc,0), which is impossible. 

If k is even, the solutions of have an even number of zeroes. So, they take 
values of the same sign at and near +00. If they were not invariant under 
the Kelvin transform, we would find two new branches, C^, bifurcating from 
(/ifc,0), which is again impossible. 



Step 3: Asymptotic behaviour of the branches 

Non trivial branches of radial solutions are contained in the region > 8}, 
that is > 2 (see Theorem [2] and the remark immediately afterwards for 
N = 2). Furthermore, with the notations of Theorem [21 by (I14p . there exists 
a unique a{N) G M such that a{a) > 2N for all a < a{N). For > 2, 
since + 2 < 2A^ < 2{N + 1), if a{a) = iV + 2, then a > a{N). Hence 
/(O) = u{0) - M^(0) > a{N) - u*j^{0) with u = Ua given by (fT5|). 

As a consequence, the branches are locally bounded for fi G [8, +00) for 
any > 2. By Step 2, is also locally bounded for /i G [8, +0x3). Since 
non trivial branches do not intersect, C^, k > 3, are all locally bounded for 
fi G [8, +00). □ 



As a simple consequence of Theorem [T21 we have the following corollary, which 
is the counterpart of Theorem [T] written for fl2^ . 

Corollary 1 For all k > 2, for all n > fi^ = 2k{k+l), there at least 2{k— 2) +2 
distinct radial solutions of (24\ ), one of them being the zero solution. 

The bifurcation diagram obtained for equation (see Fig. 1, left) is eas- 
ily transformed into a bifurcation diagram for the solutions of ([HD with A = 
2'k{N + 2) (see Fig. 1, right) through the transformation u = f + m^. In 
the case of equation ([S]), branches bifurcate from the set of trivial solutions 
C ■= {{N, \ log(2(A^ + 2))), in the representation (A^, a = u{0)). 

Based on numerical evidence (see Fig. 1), it is reasonable to conjecture that, 
in contrast, the branch C^=2 admits a vertical asymptote in the sense that 
as s +00, then A^ converges to 2, which is the only admissible value by 
f fMl) . So for (/i(s),/s) G C2, fs should develop a concentration phenomenon 
at the origin, and as s ^ +00, we should have: a +00, A^ — 2+ and 
(-i+\ %\2 G^^" Stt 5z=q, weakly in the sense of measures. 
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Fig. 1. Bifurcation diagram in the representation {N,f{0)) for equation \24^ fl^ft) 
and {N, a) for equation ^ with X = 27r(A^+2) (right). Non trivial branches bifurcate 
from Nk = 4, 10, 18, 28,. . . 

3.4 Non radially symmetric solutions 

As already described in Section 12.11 to any solution u 7^ u^j of ( ITTl) such that 
/ig2 (1 + e^" dx = 2n {N + 2), we can associate a punctured sphere of non 

radially symmetric solutions of (fTTl) . Me with e G §^ \ {N,S}, satisfying also 
1^2 (1 + Ixp)^ e2«= dx = 2tt{N + 2) for all e G \ {N, S}. And so, for > A^^, 
there are at least 2{k — 2) + 1 punctured spheres of non radially symmetric 
solutions to dm) at level A = 27r(A^ + 2). 



4 Further results, numerical observations and conjectures 

In the study of ( |T5l) . multiplicity results for general values of a are difficult 
to deduce from Proposition [3] since they require a detailed analysis of the 
nature of each critical point: maximum, minimum, and even more in the case 
of an inflection point, as well as precise estimates of the corresponding critical 
value. Numerically, Proposition [3] gives straightforward results, which can be 
observed directly from the plots of the curve a \—>- a{a) for various values of 
A^, see Figs. 2, 3. 

Let us give some details. We consider the solution Ua of (fT2l) parametrized 
by a = Ma(0). Recall that for a given A^ > 1, the curve a 1— >• a{a) = /o°°(l + 
^2-^N^2ua jg g^p]^ lima^_ooa(a) = 2(A^ + 1), lima_,ooa(a) = 2A^, and 
for A^ large enough, its range is an interval [a at, 2(A^+ 1)) if is achieved, or 
(2A^, 2(A^ + 1)) otherwise. If A^ > 2, ajv := min„6Ka(a) < a{a*^) = N + 2 < 
2N, where := ^ log(2 (A^ + 2)), which provides a multiplicity result for ( |T5i) 
in the range a e {a^, 2N) D {N + 2, 2N). 

As a function of A^ > 0, we observe that = infag^ a (a) < 2A^ if and only if 
A^ > A'o, where A"o is numerically found of the order of 1.27 ± 0.02, although 
its exact value is not easy to determine, see Fig. 4, left. In the range (A'o, 20), 
we observe that ttAf is achieved by the first critical point, see Fig. 4, right. 
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Fig. 2. Curves a a{a) for various values of N : N = 25 (left) and N = 1, 2, 
3,. . . 12 (right). The point (a^, + 2) corresponding to the explicit solution [T6\) is 
represented by a gray dot. 




Fig. 3. Curves a i— > a(a) for N = N^, A; = 2, 3, 4, 5 (left) and a i— > a{a) — 2N for 
N = 1, 3, 5,. . . 19. The function N i-^ on — 2N is monotone decreasing. 




Fig. 4. Critical points (left) and critical values (divided by AN, right) of a, as a 
function of N . 

For any N > NQ,we observe that there are at least two radial solutions of (ITS!) 
for any a G {a^, 2N). By the observations of Section it follows that there 
should be at least a radial solution and a whole punctured sphere of non 
radially symmetric solutions at level a = + 2. This supports a conjecture 
by C.-S. Lin in [25J. Actually, we can state the following result, which slightly 
improves on Theorem [2], (iii), and rigorously defines No. 

Proposition 13 There exists Nq G (1,2) such that, for all N G {Nq,oo), 
ajsf = miUagiR Q;(a) < 2A^, and for all a G {a^, 2N), there exists at least two 
solutions of /[T5\) . 

Proof. The proof relies on the continuity of the curve a i— > Q;(a) with respect 



19 



to and the fact that ^'(ag) < 0, thus proving that < 2A^ also for N < 2, 
close enough to 2. □ 



It seems that as ^ A^^q, > A'q, a at is achieved by a unique oo, so 

we may conjecture that / fl5j) has multiple solutions for a G (aAr,2A^) if and 
only if N G (A^o,oo), for some Nq G (1,2); for N G (0,A^o), the solution is 
unique, whenever it exists. A possible way to tackle such a conjecture could 
be to show that for A^ > 1, the function A^ i— > — 2N is monotonically 
nonincreasing in A^, as it appears to be the case in our numerical study, see 
Fig. 4, right, and to exploit the fact that oat — > cx) as A^ — A'o, A^ > A'o, see 
Fig. 4, left. 

On the basis of our numerical results, we may also conjecture that for A'o < 
A^ < 10 and a G [a^, 2A^), there exist exactly two radially symmetric solutions 
of ([8]). This conjecture is supported by the bifurcation analysis of Section [3l3] 
concerning the specific value Q; = A^ + 2G(a7v, 2A^) for N > 2 and A^ 7^ 4. 
Note that for A^ = 4, ajv = A^ + 2 should hold. As A^ increases, the curves 
a a{a) appear to have more and more critical points. Thus, for suitable 
values of a, the number of solutions increases as A^ increases. We have already 
checked this fact in Theorem [1] for X = 2n a, a = N + 2, but apparently it 
also holds for other values of a. 

The last observations are concerned with the function a \—>- Jn{o) defined by 
f l22l) . It seems that such a function always takes positive values on {—00, c{N)) 
and negative values on (c(A^), 00), for some c{N) > 0. See Fig. 5, left. We may 
formulate this as the following conjecture: There exists a function N ^-^ c{N) 
on (0, +00) such that Jn{o) = if and only if a = c{N) and JAr(a) > if and 
only if a < c{N). See Fig. 5, right. 




10 20 30 40 50 



Fig. 5. The function a 1-^ </Ar(a) for N = I, 3, 5, ...11 (left) and the curve 
N I— > c{N), where, at N fixed, c{N) is the first positive zero of a ^ Jiy{a); the 
dotted line corresponds to N ^ Iog(2(A^ + 2)). These two curves are tangent at 
N = W = N3 and N = 28 = N5 (right). 

Quite interesting is the comparison of with c{N). By Proposition [TOl we 
know that j{k) := JN{ci*Nj = if A; is odd and j{k) > if A; is even. Recall that 
a'{a*j^) = if and only if A^ = A^^. We observe numerically, and conjecture, 
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that: For any N > 2, a*j^ < c{N), with equality if and only if N = N21+1 for 
some I > 1. This is observed numerically with a very high accuracy for A; = 3, 
5, 7, 9, 11, see Figs. 5, right, and also 6, left. 

Summarizing the results of Section 13.21 we have shown that, in the interval 
(c(iV), 00) 3 a, the number of nodes of (pa given by ( JTTl) increases as grows 
each time a new critical point of a appears. This needs to be interpreted in 
terms of Morse index, which is still an open question. 

To investigate whether a critical point of a is a local minimum, we may look 
at the functional 

+00 

KN{a) := J {l + rYe^''^{iJa + 2ipl)rdr 


where ipa solves the ordinary differential equation 

' < + ^ + 2 (1 + rYe^""^ {ilJa + 2^l) = 0, re (0, +00) , 

< 

_Mo) = o, ^1(0) = 0. 

We have indeed a"{a) = 2Kjs[{a). No simple criterion for the positivity of 
K]\f{a) is known, but our numerical results at level a = + 2, see Fig. 6, 
right, combine very well with the results of Theorem [1] and the bifurcation 
diagrams shown in Fig. 1. 




Fig. 6. Left: the curve N 1— > JAr(a|^) is nonnegative and achieves its minimum 
value, 0, (resp. local maxima) for N = N21+1, I > 1 (resp. N = N2i). Right: the 
curve N 1-^ Ki\i{a*j^) changes sign, but is always nonnegative when a'{N) = 0. When 
^ = ^21, I ^ 1; Kiy{a*j^) is positive. 
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